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LONG-TERM  GOALS 

This  study  contributes  to  our  long-term  efforts  toward  understanding: 

•  Mixed  layer  dynamics 

•  Processes  that  communicate  atmospheric  forcing  to  the  ocean  interior 

•  Frontal  dynamics,  in  particular  the  role  of  surface  forcing  in  lateral  mixing 

•  The  interaction  of  finescale  and  submesoscale  upper-ocean  mixing  at  fronts. 

OBJECTIVES 

Physically-based  parameterizations  of  vertical  mixed  layer  fluxes  in  ocean  models  characterize 
turbulent  mixing  at  length  scales  smaller  than  the  layer  depth,  but  neglect  the  dynamics  of  unresolved 
horizontal  mixing  processes  below  their  O(l)-O(10)  km  horizontal  resolution  scale.  Numerical 
modeling  can  be  used  to  test  new  and  existing  scaling  predictions  of  surface  boundary  layer  horizontal 
mixing  processes  in  regions  of  significant  horizontal  variability,  as  commonly  found  in  major  ocean 
fronts  and  coastal  regions.  The  goal  of  numerical  modeling  work  in  in  this  year  of  the  DRI  is  to 
quantify  relationships  between  surface  fluxes  of  heat,  energy  and  momentum,  the  available  baroclinic 
potential  energy,  the  resultant  vertical  mixing  and  geostrophic  imbalance,  and  the  ensuing  dependence 
of  lateral  mixing  at  successively  larger  scales  on  atmospheric  forcing.  The  objective  in  simulation 
analysis  and  model-data  comparison  is  to  develop  dynamic  scalings  for  the  horizontal  and  vertical 
components  of  turbulent  kinetic  energy  and  fluxes  in  baroclinic  upper  ocean  environments.  These 
scalings  will  be  used  to  compare  with  observations  of  lateral  mixing  and  to  tune  modifications  of 
existing  upper  ocean  boundary  layer  models. 

APPROACH 

Large  Eddy  Simulation  techniques  are  employed  in  the  context  of  a  background  lateral  density  gradient 
where  perturbations  are  doubly  periodic.  This  baroclinic  column  geometry  is  used  to  simulate  the 
response  of  turbulence  on  scales  smaller  that  the  domain,  and  simulations  of  3D  large -eddy  turbulence 
in  boundary  layers  of  depth  10m<  7/ML<100m  can  be  compared  with  in  situ  measurements  of 
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turbulence  and  dispersion.  Model  forcing  cases  encompass  mean  buoyancy  gradients  4x10'  s' 
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<M  <2x10'  s'  ,  destabilizing  surface  buoyancy  losses  to  100  W/m  ,  and  winds  5m/s<Uio<18m/s  both 
with  and  without  wave  forcing  through  the  Craik-Leibovich  mechanism.  Domains  are  1-3  km  wide 
across  fronts,  0.4-1  km  across-front,  and  50-100m  deep  to  accommodate  mixed  layers  40<H<80m,  and 
isotropically  resolved  to  1-2  m  using  viscosity  based  on  an  advected  subgrid  TKE  budget.  Simulations 
can  incorporate  virtual  Lagrangian  Floats,  gliders  and  drifters,  to  provide  a  basis  for  interpreting  these 
small-scale  mixing  measurements. 

WORK  COMPLETED 

Numerical  modeling  in  this  year  has  been  motivated  by  differences  between  the  ~400-600m  scales  of 
coherent  lateral  structures  in  prior  simulations  of  baroclinic  mixed  layer  turbulence  in  a  front  of  the 
California  current  observed  during  the  AESOP  DRl  (Fig.  1),  and  the  \-3km  scales  predicted  for  the 
lateral  length  scales  in  wind-driven  turbulence  with  symmetric  instabilities  (SI)  by  in  results  for  similar 
LES  cases  of  Thomas  and  Taylor  (2010)  and  of  Taylor  and  Ferrari  (2010).  Numerical  simulations 
carried  out  this  year  include  forcing  and  domain  cases  to  reproduce  those  reported  in  these  studies,  and 
additional  cases  at  higher  and  lower  wind  stress  and  surface  buoyancy  loss  forcing,  as  well  as  at  larger 
and  smaller  imposed  background  baroclinic  density  gradients.  Ongoing  simulations  are  testing  the 
effect  of  the  short  down-front  dimensions  of  the  oblong  domain  geometry  on  realized  submesoscale 
instabilities. 


Surface  Buoyancy  (m/s^) 


400  600  800 

X-front  distance  -  East  (m) 


Figure  1:  Surface  buoyancy  b  from  LES  of  baroclinic  mixed  layer  in  the  California  current  with 
mean  lateral  stratification  M^=6.3xl0'^s'^,  x-y  periodic  perturbations,  and  time-dependent  wind, 

waves  and  heat  fluxes. 
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RESULTS 


Results  previously  reported  for  baroelinic  mixed  layer  eases  forced  by  steady  surface  buoyancy  loss  or 
down-front  wind  stress  were  generally  reproduced,  but  with  significant  quantitative  differences. 
Cross-frontal  length  scales  predicted  to  grow  under  the  effect  of  SI  to  --Skm  in  wind-driven  baroclinc 
layers  with  M  =6.3e'  s'  and  surface  stress  0. 1  N/m  grew  only  to  ~lkm.  For  SI  driven  by  buoyancy 
flux  the  slope  of  SI  was  reproduced,  but  with  cross-frontal  wavenumbers  2-3  times  higher. 


a)  Wind  only:  crossfrontal  section  of  crossfront  U  (cm/s) 
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b)  Wind  only:  downfront/y-  average  of  crossfront  anomaly  U'  (cm/s) 


d: 


c)  Wind  only:  Buoyancy  averaged  downfront  ( 6  b  =  0.05  cm/fe^) 
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Figure  2:  LES  of  baroelinic  front  forced  by  wind-only:  a)  Cross-frontal  x-z  section  of  cross-frontal 
velocity  U;  b)  downfront  (y-averaged)  cross-frontal  velocity  anomaly  <U’>y  ;  c)  contours  of 
isopycnals  and  mean  buoyancy  and  <b-bo>y  relative  to  central  surface. 


The  slow  growth  of  lateral  length  scales  due  to  SI  appears  to  account  for  the  remaining  differences 
between  steady-state  simulations  and  unsteady  forcing  cases.  Figure  2  illustrates  the  cross-frontal 
velocity  structures  sloping  along  isopycnals,  after  a  day  of  steady  forcing,  with  ~500m  lateral  scales, 
comparable  to  those  found  for  the  California  current  with  significant  diurnal  and  higher  frequency 
variability  (Fig.  1). 
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time  (Inertial  scale,  2;i/f ) 


Figure  3:  Evolution  of  front  scales  in  cross-frontal  surface  stratification  Mq^  over  the  course  of 
simulations  forced  by  (a)  0.1  N/m^  wind  stress  and  (b)  with  additional  forcing  from  a  fully  developed 

wave  field. 

The  evolution  of  surface  front  structures  mirrors  that  of  cross-frontal  circulations  due  to  SI  in  the  lower 
layer  for  simulation  both  with  and  without  surface  wave  forcing.  While  the  widths  of  individual  fronts 
is  altered  by  the  Craik-Liebovich  vortex  force,  the  growth  of  inter-frontal  scales  (Fig.  3,  4a)  and  lower 
layer  stratification  remains  very  similar  between  the  two  cases.  The  growth  to  0(1)  km  scales  of  fully 
developed  SI  after  3  inertial  timescales  in  both  cases  corresponds  to  lower  the  growth  and  equilibrium 
bulk  Richardson  numbers  between  0.3-0. 4  (Fig.  4b). 
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b)  Bulk  Richardson  Number 


c)  Horizontal  u  ,w  &  Diss=  [  1  rAz 
'  rms’  rms  '  ‘ 


time  (Inertial  scale,  2x/f ) 


Figure  4:  Evolution  to  steady-state  of  (a)  inter-frontal  length  scale,  (b)  bulk  Richardson  number  and 
(c)  u* -scaling  of  horizontal  (UrmJ  and  vertical  (w^J  TKE  components  and  vertically-integrated 
dissipation  (‘Diss’),  for  LES  cases  with  and  without  surface  waves. 


While  mean  dissipation  (Fig  4e)  is  elevated  several-fold  by  the  additional  TKE  production  from 
redistributing  surface  momentum  flux  down  the  geostrophic  shear  in  the  mixed  layer,  also  known  as 
the  ‘Ekman  buoyancy  flux’.  Mean  vertical  TKE  Wrms  is  not  correspondingly  increased,  except  (Fig.  5) 
locally  in  rough  correspondence  to  the  variability  in  surface  M^,  with  a  similar  pattern  for  dissipation. 
Horizontal  TKE  components  Urms^,  Vrms^  increase  steadily  to  full  SI  levels. 
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Wind  +  Waves:  Vertical  turbulence  scale  w  /u* 
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Eigure  5:  Top:  Section  of  vertical  TKE  scale  Wrms^  from  along-front  averaging,  fully  developed  SI 
state.  Bottom:  Surface  cross-frontal  buoyancy  gradient  Mo  scaled  on  its  average  value 
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The  results  from  these  and  other  completed  simulations,  and  from  those  currently  under  way,  are  being 
analyzed  to  determine  the  scaling  behaviors  of  lateral  kinetic  energy  and  flux  components.  While  the 
behavior  of  LES-domain-averaged  dissipation  and  vertical  kinetic  energy  is  not  consistent  with 
observations  of  Sl-driven  baroclinic  boundary  layer  turbulence  observed  in  the  Kuroshio  (D’ Asaro  et 
al,  2011)  and  in  the  Gulf  Stream  during  LatMix  2012  field  experiments,  turbulence  scaling  based  on 
conditional  averaging  at  fronts  in  these  LES  results  do  appear  consistent  with  the  observations  for 
downfront  wind  cases,  provided  that  scaling  front  interior  dissipation  and  Wrms  with  ‘Ekman  Buoyancy 
Flux’  (EBF)  forcing  (Thomas,  2005)  is  not  additive  with  forcing  by  wind  stress  and  surface  waves 
(Fig.  6).  A  case  with  wind  and  waves  directed  partially  cross-front  is  notably  not  consistent  with  EBF 
scaling  expectations. 


Scaling  dissipation  with  EBF,  domain  averages  (stars)  and  conditional  averages  (dots)  in  fronts 
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Figure  6:  Scaling  with  EBF  of  layer-integrated  dissipation  <s>H  on  surface  friction  velocity  u*, 
domain-averaged  (stars)  and  averaged  conditionally  on  bins  of  surface  horizontal  buoyancy 
gradient  (dots).  EBF  is  the  dot  product  of  Ekman  transport  and  surface  buoyancy  gradient.  LES 
case  with  wind  &  waves  directed  30  degrees  toward  the  lighter  side  of  front  does  not  conform  to  1:1 

scaling  line  with  EBE. 


IMPACT/APPLICATIONS 

Eateral  Mixing  DRI  results  bear  on  the  predictive  skills  of  regional  scale  models  with  0(1-10)  km 
resolution.  At  these  scales  the  parameterizations  of  both  vertical  and  lateral  fluxes  are  not  well 
understood  or  tested,  especially  in  energetic  frontal  environments,  or  in  subsurface  environment  where 
lateral  mixing  is  driven  by  the  relaxation  of  stratification  anomalies  produced  by  turbulent  mixing 
events. 
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